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Abstract 

We optimize the turning on of a one-dimensional optical potential, Vl{x, t) = 5(t)Vo cos^(A;x) to 
obtain the optimal turn-on function S{t) so as to load a Bose-Einstein condensate into the ground 
state of the optical lattice of depth Vq. Specifically, we minimize interband excitations at the end of 
the turn-on of the optical potential at the final ramp time tj., where S{tj.) = 1, given that ^(O) = 0. 
Detailed numerical calculations confirm that a simple unit cell model is an excellent approximation 
when the turn-on time t^ is long compared with the inverse of the band excitation frequency and 
short in comparison with nonlinear time h/ where ^ is the chemical potential of the condensate. 
We demonstrate using the Gross-Pitaevskii equation with an optimal turn-on function S{t) that 
the ground state of the optical lattice can be loaded with very little excitation even for times tr on 
the order of the inverse band excitation frequency. 

PACS numbers: PACS Numbers: 3.75.Lm, 03.75.-b, 67.90. +z, 32.80.Qk 
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I. INTRODUCTION 



A number of competing schemes for implementing quantum information and quantum 
computing are being explored; five physical systems have been proposed as quantum logic 
gates: ion traps high-Q optical cavities P, 0], nuclear magnetic resonance systems j^, 
solid-state qubits (semiconductor quantum-dot and Josephson-junction devices) P, and 
ultra-cold neutral atoms in optical lattices [3, The proposals for using atoms in optical 
lattices can be implemented by first loading Bose condensed atoms that are held by a weak 
magnetic field into an optical lattice by gradually turning on the optical potential to its 
desired strength. Upon increasing the intensity of the optical lattice to a critical intensity, the 
Bose-Einstein condensate (BEC) will undergo a quantum phase transition from a superfiuid 



state to a Mott-insulator state P]. One thereby can obtain one atom per lattice site in the 
ground state of the system; these atoms can serve as qubits. This suggestion of preparing a 
Mott-insulator state has recently led to a seminal experiment wherein the quantum phase 
transition was observed. In principle, starting with a BEG in a trap and turning on an optical 
lattice of sufficient well depth in a sufficiently adiabatic manner, prepares the Mott-insulator 
state. In practice, it is easy to turn on the optical lattice adiabatically with respect to band 
excitation (excitation from one band to another); however, it is substantially more difficult 
to turn on the optical lattice adiabatically with respect to quasi-momentum excitation. The 
second, more stringent form of adiabaticity requires that the optical lattice be switched 
on slowly with respect to mean-field interactions and tunneling dynamics between optical 
lattice sites, and hence typically requires milliseconds Q]. We refer to the first form of 
adiabaticity as 'mterband adiabaticity' and the second form as 'miraband adiabaticity'. 
The intraband adiabaticity condition has been demonstrated in one-dimensional lattices by 
Orzel et al. and ultimately led to the pioneering experimental demonstration of the 
Mott-insulator transition [l2|. 

The goal of this paper is to investigate the loading of atoms into an optical lattice in as 
short a time as possible, while maintaining interband adiabaticity, so as to obtain maximal 
atomic population in the ground state band with minimal band excitation. We consider 
experiments of the type described in Ref. 3]; which measured the interband nonadiabaticity. 
We show how to optimize the turning on of a one-dimensional (ID) optical lattice so as 
to minimize the interband nonadiabaticity. We carry out calculations using the Gross- 
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Pitaevskii (GP) equation, as well as simplified models. We show that an optimized lattice 
turn-on results in very low nonadiabaticity even for loading times tr comparable to l/Acj, 
where hAu is the excitation energy of the first band that can be excited. On such timescales, 
the nonlinear term in the GP equation is small and has little effect on the dynamics. 

Section II discusses the theoretical models used, and describes how to determine the 
optimized turn-on of the lattice. It also describes a very simple unit cell model that is in 
excellent agreement with the full calculations. Section III describes our results and compares 
to experimental data. A final Section summarizes our conclusions. 



II. THEORY 



Our theoretical investigation was stimulated by experimental work to load a BEG into 
a one-dimensional optical lattice so as to create the ground state of atoms in the optical 
lattice 3|- We first briefly describe the optical potential and its turn-on, then consider the 
mean-field description of the process, develop the simple unit-cell model, and finally describe 
the optimization procedure used. 



A. The Optical Potential 

We start with a condensate of N atoms in a trap with frequencies u^, Uy, and u^. At time 
t = 0, a ID optical lattice with potential 

Vl{x, t) = -S{t)Vo cos2(A;x) (1) 

is turned on along the x direction. The lattice periodicity equals ir/k, where k = ujl/c and 
ujl/2ti are the wave vector and frequency of the laser beams that form the optical lattice. 
The lattice beams are sufficiently far-detuned so that spontaneous emission is negligible on 
the time scale of the dynamics described here. The parameter Vq is the final lattice depth 
after completion of the turn on process at ramp time tr- The ramp turn-on function S{t) 
has only two constrains: S'(O) = and S{tr) = 1. For diagnostic purposes the lattice is held 
on with fixed depth Vq until time tf = tr + thoid, then turned off to allow for condensate 
expansion and imaging. Fig. shows two examples of ramp functions with tr = 20 /xs: a 
piece-wise cubic function with function and derivative matched at tm = 16 /is and a fit to 



the experimental ramp function used in Ref. 



3- 




B. Mean-field Description 

The GP equation for the BEG wave function \I/(r, t) (the order parameter of the conden- 
sate) is 

ih— = v^^ + v(v,t)-^ + m^m , (2) 

ot 2m m 
where V{T,t) = Vi;rap(r) for times t < and V{r,t) = VL{x,t) for t > 0, and the nonlinear 
term is proportional to the two-body s-wave scattering length a, and m is the atomic mass. 
We assume the trap is turned off at t = 0. The lattice turn on time tr is assumed to be very 
short compared to the time needed for the condensate to expand significantly. 

The initial condensate wave function at t = is that for a trap with time-invariant poten- 
tial Vtrap(r)- Although the momentum distribution of the initial condensate wave function 
\&(r, 0) is sharply peaked close to zero, the final condensate wave function \&(r, t^) in the 
lattice has developed momentum components peaked near 0, ±2hk, ±Ahk, .... These com- 
ponents appear even if only the ground state band of the lattice is occupied. Gonsequently, 
the subsequent free evolution of the condensate after tf results in the further physical split- 
ting of the condensate wave packet into spatially separate parts having central momenta 
2nhk, where n = 0, ±1, . . . . The experiment consists of measuring the fraction Pnk{tf) of 
the initial total number of condensate atoms that appear in the nth wave packet. As we 
describe below, the signature of nonadiabatic excitation during the turn-on of the ramp is 
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oscillatory behavior of Pnk{tf) with varying hold times thoid- 

The central goal of this paper is to determine the optimized ramp function Sopt{t) for 
a given tr and Vq, i.e., the function S{t) which leads to minimal oscillatory amplitude of 
Pnk{tf), and hence minimal excitation of the final condensate wave packet at time tr (see 
Eq. ©). Our optimization will be carried out using the GP equation. We do not require 
that our system remain adiabatic during the whole time interval between t = and tr, but 
only that the inter-band excitation be as small as possible at the final time tr- It should be 
noted that the mean-field treatment used here can not be used upon increasing the optical 
lattice potential to the point where the BEG undergoes a quantum phase transition from 
its BEG-like superfluid state to a Mott-insulator state. Once the number of atoms in each 
optical potential well becomes small (i.e., no longer large compared with unity), a field 
theory description must be used. 

It is easy to calculate the experimental observables from the GP wave function \l/(x, tj). 
By performing a Fourier transform it is easy to calculate the fraction Pnk{tf) of atoms 
associated with each sharply peaked momentum component n. The degree of interband 
non-adiabaticity /nonad is determined experimentally by the amplitude of the oscillations in 
Pok{tf) versus tj: 

/nonad ~ POk,msix POfc,min ; (3) 

where po/c,max and pofc,min are the maximum and minimum values of Pnkitf) for n = over an 
interval tf — tr = thoid which is chosen to be long enough to include at least one oscillation 
of Pok{tf)- If there is no interband excitation, then /nonad = 0. 

Figure |21 shows various components of the BEG energy as they evolve in time when an 
optical potential with 2Tx/k = 590 nm, and Vq = 14£'r, where the recoil energy En is defined 
as Er = h^k'^/2m, is turned on linearly over 20 /is, held constant for 20 /iS, and finally turned 
off over 20 fis. The trap with frequencies = 84 Hz, Uy = 59.4 Hz, and z/^ = 42 Hz holds 
10^ Na F = 1,M = -1 atoms, for which a = 2.8 nm and Er/H = 25 kHz. Notice that 
the total energy of the system is constant when the potentials do not depend on time. The 
mean field energy remains negligible in comparison to the kinetic and potential energies 
that result from turning on the optical potential. Gonsequently, calculations with 10,000 
and 1x10^ atoms are indistinguishable on the scale of Fig. El We also checked that on these 
short timescales the calculated nonadiabaticities from both 3D and ID versions of the GP 
equations are indistinguishable. Thus, ID linear Schrodinger dynamics should provide an 
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FIG. 2: Kinetic, optical, mean-field and total energies versus time starting from a condensate with 
10'* Na atoms in a trap with frequencies of 84 Hz, 59.4 Hz, and 42 Hz in the respective x, y, and 
z directions. A 20 /US linear ramp is applied ending at a latttice depth of Vq = 14£^/j followed by a 
20 //s hold time and a 20 //s linear ramp down. 

excellent approximation to the nonadiabatic dynamics. 
C. Unit-cell Model 

We only consider relatively short time scales t,. of lattice loading so that mean-field 
interactions do not play an important role, as shown in Section FlI Bl We assume t,. is short 
compared with the nonlinear time ^nl = ''^/(4vr/ia|\l/peakP) = ^//-tj where |\l/peakP is the 
peak value of |\E'p at t = 0, and yU is the chemical potential. Consequently, we can develop 
a simple linear 'unit cell' model of the loading process that is in excellent agreement with 
the full GP dynamics. In the unit cell model, the natural units of length, energy, and time 
are l/Zc, the recoil energy Er = and a characteristic lattice band excitation time 

^band = 271 / Au respectively. 

Fig. El shows the Wigner-Seitz primitive unit cell of an infinite optical lattice of depth 
Vo = 14i?/j. Solving the Schrodinger equation for the potential —Vq cos^ kx with periodic 
boundary conditions on the interval —Ti/2k < x < Ti/2k defines a series of eigenvalues 
EiiVo) and eigenfunctions |z, Vq) that describe the unit cell wave function for the g = edge 
of each band i, where q is the lattice momentum Q]. The boxes in Fig. El indicate the width 
of the energy bands. The width of the lowest band is too narrow to be observed on the 
energy scale of the figure. The width of the third band extends beyond the top of the lattice 
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FIG. 3: Unit cell lattice potential Vl{x) versus x, showing the energy spread of the first three 
energy bands for a Na lattice with 27r/A; = 590 nm and Vq = 14£'r, where Efi/h = 25 kHz. 
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FIG. 4: Energies of the first three q = levels versus lattice depth Vo/h for the unit cell shown 
in Fig. 121 . The second and third levels are degenerate with an energy 4£'/j in free space, but 
split into symmetric {E^) and antisymmetric {E2) states with increasing lattice depth. Turn- 
ing on the lattice rapidly couples the symmetric i = 1 ground and i = 3 levels, but not the i = 2 one. 

potential. Fig. |3] shows the energies Ei for the q = levels for the first three lattice bands 
as the lattice depth Vq is increased. Since symmetry considerations only allow excitation 
of odd numbered bands from the first n = 1 band during lattice turn on, the relevant trap 
excitation parameter Auj is [E-^ — Ei)/h. For the example in Fig. |2l Auj = 27r(250kHz) and 

tband = 4 /iS. 

An adiabatic basis set expansion can be used to calculate the time-dependent amplitudes 
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aj{t) of the instantaneous eigenvectors during the dynamics. The basis set expansion is 
taken in the form, 

vl/(t) = 5^a,(t)e-/o^.(*')'^*'/^|j(t)), (4) 
j 

where \j(t)), j = 1,2,..., is the instantaneous eigenvector with energy eigenvalue Ej(t) 
corresponding to an optical potential strength S{t)Vo. The nonadiabatic dynamics is then 
calculated from the coupled set of equations, 

^ ^E,{t)-E,{t) 

where we have defined the coefficient aj{t) = — Ej{t')dt' /h to simplify the notation. 
The even index and odd index eigenstates are not coupled because the potential and its 
time derivative are an even function of x and the odd (even) eigenstates are symmetric 
(antisymmetric). We numerically calculate Ei{t), Ej(t) and {i(t)\H{t)\j(t)) , and solve the 
coupled differential equations by using a variable step size integrator. We will call this 
method the converged unit cell model. 

We have verified by numerical calculations that a truncated basis set expansion with the 
lowest two terms is a good approximation unless the ramp time t,. becomes comparable to 
^band or less: 

vl>(t) = ai(t)e^"i«|l(t)) + a3(t)e'"^(*)|3(t)), (5) 
This leads to the following simplified 2-equation unit cell model: 

An even simpler one-equation approximation, which is acceptable for sufficiently large tj., 
is to set ai = 1 + iO and solve for a-sit) by integrating Eq. ((Tj), so that 

J, E,{t')-E,{t'f 
We will call this the one-equation unit cell model. 
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The fraction of the total number of atoms in the wave packet with mean momentum nk, 
Pnk{t), can be calculated in terms of the amplitudes aj{t): 

2 



Pnkit) = |(2nA;|*(t))p 



J2<^At){'2nk\j{t))e 



iaj (t) 



^ |ai(t)(2nfc|l(t))e'"i(*) + as{t){2nk\3{t))e''''^'^\\ (9) 
When n > 0, the total population with magnitude of momentum \nk\ is 

P\nk\{t) = Pnk{t) +p_nk(t) . (10) 

The observable /nonad in Eq- © can be calculated from Pokit). Using the one-equation 
model of Eq. (jH)), where Oi = 1, /nonad simplifies to 

/nonad = 4(0fc|l)(0A;|3)|a3|. (11) 

The expression in Eq. fllip is evaluated for t > tr, where none of Vl, |3), or |a3| is 
changing in time. 

The fidelity F of the ramp can be defined as the fraction of population that is left in the 
ground state after the ramp: F = |aip, where ai(t) is evaluated for a time after tr- In the 
two-state model, F = 1 — |a3p. Using the result in Eq. (jllj) . we can get a relation between 
F and /nonad: 

^ " ^ ~ 16|(0A;|1)(0A;|3)|2^"°""^ ^^^^ 
For the case of the Vq = I^Er Na lattice discussed in the figures, F = 1 — 0.28|/nonadP- 

Fig. El compares the numerical solutions for the two-equation model of Eqs. © - © and 
the one-equation model of Eq. © for a Vq = 14£'/j lattice with the cubic ramp function in 
Fig. ^ Clearly, the simpler one-equation model provides a very good approximation in this 
case. Figure El shows that the degree of nonadiabaticity at the end of the ramp can be much 
smaller than during the turn-on portion of the ramp. This suggests that we might be able 
to optimize the ramp so that the final value of jasp at the end of the ramp is small, even if 
it becomes large during the ramp. 

Figure El shows the populations of the free space momentum states, pok(t) and p\2k\(t), 
for several different calculations. The Vq = 14i?K lattice is turned on with the cubic ramp 
function in Fig. ^ The figure shows the momentum state populations that correspond to 
the adiabatic eigenstates, |(2n/c|j)p, as well as those obtained with dynamical calculations 
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FIG. 5: Time evolution of the real and imaginary parts of 03 calculated with the 14ii^/j deep, 20 
/is cubic ramp shown in Fig. ^ The results are shown for the two equation model and the one 
equation model. Eqs. (jH)-© and Eq. Q respectively. 
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FIG. 6: Population of the Ok and \2k\ momentum components versus time. The optical lattice is 
turned on to a depth of 14£'/j with the 20 /is cubic ramp shown in Fig. ^ 

based on the GP equation and the two-equation unit cell model of Eqs. © - ©• The figure 
demonstrates the very good agreement between the GP and unit cell model. The oscillations 
of the dynamical populations around the adiabatic ones are due to nonadiabatic excitation 
due to the fast ramp. The amplitude of these oscillations correspond to /nonad = 0.07 and a 
fidelity of F = 0.9985. 




— Re 2 eqs. 

— Im 2 eqs. 
-- Re a^ 1 eq. 
-- Im a^ 1 eq. 
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D. Optimization 

The oscillations evident in Fig. IHl are due to the small but finite value of nonadiabatic 
excitation that reduces the fidelity F of the ramp. We will use an unconstrained nonlinear 
optimization to minimize the degree of nonadiabatic excitation at the end of the ramp and 
thus maximize the fidelity F. For this purpose we use an A^th order polynomial form for the 
ramp function S{t), with S'(O) = and S{tr) = 1. The form of S{t) is given by 

N+l 

S{t) = c,f , (13) 

where ctv+i = (1 — '^f=QCitl) /t^~^^ insures that S{tr) = 1. We then can use a standard 
computer code to optimize the polynomial coefficients q, i = 1, . . . , N to give the minimal 
/nonad • Wc cau Optimize /nonad Calculated from the GP equation, or from the one- or two- 
equation unit cell models. We will use the GP equation, since this allows us to work with 
very short ramps, where the nonadiabaticity is large during the ramp. We will use the 
simpler models for interpreting the results. 



III. RESULTS 



A. Comparison to experiment 

In the experiment of Ref. a BEG with 3 x 10^ sodium atoms in the 3Si/2, = 1, 
rriF = —I state with no discernible thermal component is prepared and held in a magnetic 
time-orbiting potential (TOP) trap with trapping frequencies in the x, y and z directions of 
27 Hz, 19 Hz and 13.5 Hz. After production of the condensate, an optical lattice consisting of 
two counter-propagating laser beams along the x-direction is turned on. The lattice beams 
are detuned about 60 GHz to the blue of the sodium D2 line at 589 nm. The spontaneous 
emission rate is negligible on the timescale of the experiments. The polarization is linear 
and parallel to the rotation axis of the TOP trap bias field (the y-axis). The condensate is 
located in the focus of the beams, which have a 1/e^ beam diameter of about 600 fim. 

In order to directly measure the adiabaticity, i.e., the efficiency of transfer into the lattice 
ground state, an experiment was carried out in which the intensity of the optical potential 
was ramped up to a stationary lattice potential strength Vq ~ 14:Er The BEG was held 
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FIG. 7: Comparison of theory and NIST experimental data (closed circles) on the population 
of the Ok momentum component versus time for a 20 /is ramp of a Vq = 14£^/j deep optical lattice. 
The open circles shows data shifted to have the same peak magnitude as that of the calculation. 



in the lattice for a time thoid? typically between and 10 /is, before suddenly switching off 
the light. The plane- wave decomposition of the lattice wave function at the switch-off time t 
depends on the previous history of the ramp prior to time t. The populations of the various 
momentum components at t is measured by allowing the atomic cloud to expand after the 
switch-off, since each 2nk momentum component eventually separates from the others in an 
individual cloud that can be imaged. Figure shows calculated and observed oscillations 
for the population of the Ok component following a ramp of 20 /is. This beating signal has 
a small amplitude, indicating that most of the population was in the lowest band. However, 
if only the ground state were populated, there would have been no beating at all. 

Figure [7| also shows our calculation of pokit) using the GP equation with the experimental 
latticepotential having the ramp shown in Fig. ^ The figure shows the original data of 
Ref. as well as data that has been shifted in phase and magnitude to coincide with 
the calculation. These adjustments are consistent with an experimental uncertainty in the 
lattice depth Vq of about 10% and the uncertainty in the ramp function. Even without 
the adjustments, it is clear that our calculation gives the right magnitude of the effect of 
nonadiabatic dynamics as measured by the oscillations in the signal after tr- 
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FIG. 8: Optimal ramp functions for 6 fis, 15 fis, 20 fis, and 25 fis ramps for a HEr lattice. 
B. Ramp Optimization Results 

Although the degree of nonadiabaticity in Fig. [7| is small, it can be made even smaller 
by optimizing the ramp turn on function. We used a numerical optimization algorithm 
to minimize /nonad as defined in Eq. (jS)) versus the polynomial fit parameters in Eq. (fT^. 
We calculated the dynamics with the GP equation, the one-equation unit cell model, the 
two-equation unit cell model, and the converged unit cell model. The one-equation unit cell 
model works very well except for short ramps with tj. on the order of tband- Only the GP 
and the converged unit cell model work well in this limit, due to excitation of higher bands 
that are not modeled within the one- or two-equation unit cell model. 

Figure |H1 shows the optimal ramps obtained for several different ramp durations, using 

= 12 in Eq. (jl3|) . Clearly, the optimal ramp shape depends strongly on the ramp duration. 
For very short tr, the ramps are not robust in the sense that slight deviations from the 
optimal ramp shape significantly lowers the fidelity. Moreover, different optimal ramps can 
be obtained upon starting with varying initial functions for S{t). Thus, the optimization 
functional being minimized can have several local minima so that several different optimized 
functions can have similar high fidelities. 

Figure M shows a^it) versus time for the optimal 6 /is and 20 fis ramps in Fig. |S| In 
both cases, there is clearly significant excitation during the ramp, however, at the end of the 
ramps, the excitation is nearly zero. Hence, the dynamics is clearly nonadiabatic for both 
ramps, but at the final time tr, there is very little excitation, i.e., high fidelity is achieved 
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FIG. 9: Real and imaginary parts of the amplitude 03(4) versus time for 6 and 20 ^us optimal 
ramps, showing strong nonadiabatic behavior as a function of time, but ending at tr with very 
small excitation. 



for both ramps. With the optimized ramps, the fidelity obtained is much higher than that 
obtained in the NIST experiment, even with the 6 /.is ramp. However, for the optimal 6 yus 
ramp the as coefficient goes to zero at the end of the ramp with a high slope, indicating 
that the fidelity will not be robust with respect to small changes in the ramp parameters. 
Also, for small t^? the optimized ramp functions can vary dramatically and suddenly with 
time and therefore may be difficult to achieve experimentally. 

Figures ITUl and ITT] show the calculated 0/c and \1k\ populations versus time as determined 
using the full GP equation and the converged unit cell model for the optimal 20 and 6 /is 
ramps respectively. The Sif) ramp functions are also shown in the figures. The converged 
unit cell model reproduces the GP results very accurately, even for the 6 /is ramp. For the 
6 /is ramp, there is already a small 4/c population component present. Nevertheless, the 
fidelity remains high. 

IV. SUMMARY 

We calculated the dynamics of a BEG upon turning on a one-dimensional optical po- 
tential. Good agreement with experimental data has been obtained. We showed that the 
optimal ramp function S'(t), with 5(0) = and S{tr) = 1, for turning on a one-dimensional 
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FIG. 10: GP and unit cell populations versus time for optimal 20 /is ramp. The fidelity is 0.999999. 
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FIG. 11: GP and unit cell populations versus time for optimal 6 fis ramp. The fidelity is 0.9996. 

optical lattice potential, Vl(x, t) = S (t)Vo cos^ (kx) , can be chosen to minimize the inter- 
band nonadiabaticity. A simple unit cell model that is in excellent agreement with the full 
calculations has been developed for cases when the turn-on time tr is relatively long. An 
optimized lattice turn-on with very low nonadiabaticity can be chosen even for short loading 
times tr comparable to 27t/Auj, where hAu is the band excitation energy. 

Even if the loading of a BEC into the lattice without causing inter-band excitation is 
readily achievable, as shown here, unless one switches on an optical lattice very slowly, the 
optical lattice causes a spat ially varying phase to accumulate across the condensate due to 
intra-band excitation. Ref. [lO| has shown analytically and numerically that a cancellation 
of this effect is possible by appropriately adjusting the harmonic trap force- const ant of the 
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external magnetic trap that confines that BEC, thereby facihtating quick loading of an 
optical lattice for quantum computing purposes. 
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